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I. INTRODUCTION 

The description of dissipation in quantum mechanical systems is a question 
which has interested scientists since decades. Brownian motors [E 0, i, i.e. 
devices able to produce useful work out of thermal forces with the help of other 
unbiased forces, provide an ideal benchmark for the investigation of quantum 
dissipative systems Q, for two reasons. First, the interaction with a dissipa- 
tive environment plays an essential role in the performance of Brownian motors. 
Second, dissipative tunneling enriches the dynamics of quantum Brownian mo- 
tors with respect to their classical counterpart, inducing features such as current 
reversals as a function of temperature [f| . Experiments in the quantum regime 
have been reported in 0, 0] . Brownian motors belong to the family of ratchet 
systems, which exhibit directed transport in the presence of unbiased forces 
which need not necessarily have a thermal origin. 

A simple model for a quantum ratchet system is provided by the Hamiltonian 

H R =p 2 /2M + V(q) (1) 

of a quantum particle of mass M in one-dimensional space. The asymmetric 
potential of periodicity L is completely characterized by the amplitudes Vi and 
phases ipi of its harmonics in the Fourier representation 

V(q)=J2 l=1 VlCOS ( 27Tl l/ L ~Vi)- (2) 

In order to bring thermal fluctuations into the system, we let it interact with a 
bath of harmonic oscillators given by the standard Caldeira-Leggett Hamiltonian 

fl B = (l/2)V' ° Pa/m a + m a Lol (x a - c a q/m a Lol) 2 . (3) 

z — 'a—l L . 

The spectral density J(w) = (tt/2) J2a=i( c a/ m a UJ a )o~{u — u a ) fully character- 
izes the bath. We consider an Ohmic bath with viscosity r/, i.e. J(w) ~ r/uj at low 
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frequency u). An additional contribution H cxt (t) = —F(t)q to the Hamiltonian 
accounts for an unbiased driving force F(t) breaking thermal equilibrium. 

The quantity of interest is the stationary velocity — lim^oo q(t). A nonzero 
particle velocity means that work is released, although all forces acting on the 
system are unbiased. We have developed two methods described in the 

next Sections, to evaluate the stationary velocity in the model presented above. 



II. QUANTUM RATCHETS WITH FEW ENERGY BANDS 

Our first method relies on a restriction to the low-energy dynamics of the 
ratchet system. Using Bloch theorem, the ratchet Hamiltonian JT]) can be diag- 
onalized exactly, yielding energy eigenstates organized in an infinite number of 
bands. When the energies associated with the driving force and the temperature 
are much smaller than the lowest band gap, it is satisfactory to keep the few 
N-q bands lying below the potential barrier only. After rotation to the position 
eingenbasis \m,j) and use of the localized character of the retained low-energy 
states, the Hamiltonian takes the form of a multi-band tight-binding model, 



Htb=^2°° [E B £m\m,j)(m,j\+Y, " , A%™\m> , j) (m, j\ 

L<=i ( A -'™ V, j + l><m,i| + A^V, j)(m, J + 1|)] , (4) 



+ V 

* — ^m,m' — 

with on-site energies e m . The interband couplings A^™ allow for vibrational 
motion inside each potential well, whereas A^*™' f and A™*™' b allow for tunneling 
to the nearest neighboring wells. 

The dissipative dynamics in this tight-binding system, driven by a non- 
adiabatic harmonic force F(t) — Fcos(ftt), can be investigated by means of 
real-time path integral techniques. One can derive a generalized master equa- 
tion for the state populations P m j(t), in terms of the transition rates induced by 
the couplings. When the transition rates remain much smaller than the driving 
frequency f2, one can average the dynamics over one driving period and obtain 
the averaged ratchet velocity at long times, 

- t \ ^ B oo / fainter, f fainter, b\ /r-\ 

v « = L K, hm , =1 P m {r m > m -r m , m )■ (5) 

The asymptotic population p"^ — lim^oo J^'jL-oo ^m,j{t) of the band m, is also 

a simple combination of the averaged transition rates T f , and r^*^' b . 

The details of the evaluation, involving the computation of the transition rates 
up to second order in the tunneling couplings, can be found in 0]. As a result, 
the ratchet velocity reveals a nonmonotonic dependence on the amplitude and 
frequency of the driving force and on dissipation. Current inversions can be 
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obtained by changing any of these parameters. At the location of an inversion, 
the system thus experiences driving-induced localization. 

The main drawback of this approach is that the validity regime for the trun- 
cation to a finite number of energy bands excludes the interesting limit of large 
driving amplitude or temperature, related to the classical limit. This has moti- 
vated the development of another approach presented in the next Section. 

III. DUALITY RELATION FOR QUANTUM RATCHETS 

Our second method is based on a duality relation [ll[ for a quantum dissipative 
system with a sinusoidal potential tilted by a time-independent (DC) force F. 
We have generalized this relation to periodic potentials of arbitrary shape as 
given in (J2J). Here again, the evaluation of the particle velocity v-dc(F) at long 
times is based on real-time path integral techniques. One first performs an 
exact expansion of the propagator of the system in powers of the amplitudes Vi 
of the potential harmonics. After this transformation, the path integrals involve 
Gaussian integrals only, which can be evaluated. The resulting series expression, 
in powers of Vi, can be related to the series expression, in powers of the couplings, 
for the particle velocity ■Cdc(-F) in a single-band tight-binding model with non- 
nearest-neighbors tunneling couplings. This relation, valid a long times, reads 

UdcCF") = F/r] - v oc (F), (6) 

where 77 denotes the viscosity in the original system. The Hamiltonian of the 
dual tight-binding system in which vuc(F) has to be evaluated takes the form 

% = E M ,Er (A m \l + m){l\ + A* m \l)(l + m\). (7) 

A striking outcome of the calculation is the relation A; = (1/2) Vi e llpi giving the 
coupling A; to the Zth-order neighbor in terms of the amplitude Vi and phase ifi 
of the Ith harmonic of the original potential. The tight-binding system ([7]) is 
tilted by a time-independent force F as the original system. It is bilinearly cou- 
pled to a different bath of harmonic oscillators, characterized by the spectral 
density Jtb(^) = >^( w )/[l + (w/7) 2 ], thus Ohmic as J{lo) with an additional 
Drudc cutoff at the frequency 7 = -q/M associated with dissipation in the origi- 
nal system. The spatial periodicity L of the tight-binding model is related to the 
periodicity L of the original system through L = L/a, with the dimensionless 
dissipation parameter a = r/L 2 /2ttH. As a consequence, the dissipation parame- 
ter is a. = l/a in the tight-binding system. A regime of weak dissipation in one 
system is thus related to a regime of strong dissipation in the dual one. 

The particle velocity vuc(F) hi the tight-binding model can be evaluated 
from the transition rates between the tight-binding states, like in the prob- 
lem described in the previous Section. The duality relation © can then be 
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FIG. 1: Ratchet velocity (solid line, left and upper axes), and load characteristics 
(broken lines, right and lower axes) for a two-harmonics ratchet potential with ampli- 
tudes Vi = 4 V2, phases (p% — 2<pi — —it/2, and a spatial periodicity L. The load force 
is denoted by K and the amplitude of the bistable driving force by F. The vertical 
arrows show the driving amplitudes chosen for the computation of the load charac- 
teristics. The thin dotted lines are guides for the eyes. The temperature has been 
set to ksT = 0.1 AV, in terms of the barrier height AV = 2.2 Vi. The viscosity rj 
and the dynamical parameters of the system can be specified through the dissipation 
parameter a = r/L 2 /2nh and the ratio between the dissipation rate 7 = r]/M and the 
classical oscillation frequency in the untilted potential Qq — 2tt^/Vi/ML 2 . A regime of 
quantum dynamics and weak dissipation has been chosen with a — 0.2 and 7 = 0.23 ho- 



used to obtain the particle velocity vr>c(F) m the original system tilted by the 
time-independent force F. From this result, one can deduce the ratchet veloc- 
ity v-r(F) = vbc{F) + vdc{— F) in a system driven by an unbiased bistable force 
switching between the values ±F, in the limit where the switching rate is much 
slower than any other time scale of the system. Another interesting quantity is 
the velocity v^(K, F) — vuc(K + F) + vuc(K — F) in a ratchet system driven 
by a bistable force ±F and subject to a constant load force K, showing the load 
characteristic of the ratchet system. If the sign of this velocity is opposite to 
the sign of the load force, it means that the ratchet system is able to yield work 
against the load, e.g. to lift it. 

The ratchet velocity and load characteristics for a potential with two harmon- 
ics are shown in Fig. [TJ The calculation follows the lines detailed in [Io| . 
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The ratchet velocity (solid line, left and upper axes) shows several reversals as a 
function of the driving amplitude. It vanishes in the absence of driving. In this 
case the system is at equilibrium and no work can be extracted according to the 
Second Principle of Thermodynamics. The ratchet velocity also vanishes when 
the driving energy FL is much larger than the potential barrier AV", because 
the potential, which is the only source of spatial asymmetry in the system, is 
irrelevant in this regime. The figure also shows the load characteristics (broken 
lines, right and lower axes) as a function of the load force K for three different 
amplitudes of the bistable driving force. For FL = 4.70 AV (dotted line), there 
is a range of negative load forces for which the velocity is positive. The ratchet 
system is able to work against the load in this regime. A range of positive load 
forces and negative velocities is obtained for FL — 3.60 AV (dashed line). Fi- 
nally, for FL = 4.15 AV (dashed-dotted line), the load characteristic crosses the 
origin, meaning that the ratchet system is tuned into an idle regime. 

An interesting extension of this work would be the evaluation of the diffusion 
coefficient and current noise in quantum ratchet systems. 
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